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Abstract 



(j ! We discuss various aspects of vector meson production, first analysing the interplay 

between perturbative and nonperturbative aspects of the QCD calculation. Using a general 
method adapted to incorporate both perturbative and nonpertubative aspects, we show 
that nonperturbative effects are important for all experimentally available values of the 
photon virtuality Q 2 . We compare the huge amount of experimental information now 

m 

available with our theoretical results obtained using a specific nonperturbative model 
\& • without free parameters, showing that quite simple features are able to explain the data. 
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I. INTRODUCTION 



Electroproduction of vector mesons provides an interesting laboratory for studying 
the interplay between perturbative and nonperturbative QCD. Most emphasis has 
been put on the perturbative side of the calculations 0, El Q H i3- S El El llfl| ■ 
There the production process is considered to be mediated by two gluon exchange 
and the coupling of the gluons to the hadron, generally a proton, is taken from a 
phenomenologically determined gluon distribution as obtained from deep inelastic 
scattering (DIS). The coupling of the exchanged photons to the produced quarks 
is treated perturbatively. This is justified if there is a truly hard scale and the 
produced quarks stay close together. 

The emphasis put on perturbation theory is understandable given its justification in 
terms of basic elements of QCD. Nevertheless we think it is useful and even necessary 
to scrutinize also the other side of the medal, namely the nonperturbative aspects, 
and investigate its consequences. To examine the roles and magnitudes of different 
effects we have chosen an approach which starts from general expressions and inves- 
tigate the limit where perturbation theory holds. We derive and discuss particularly 
the deviations from perturbative QCD induced by nonperturbative effects. 

There are several important reasons for this approach. 

• Even if the produced vector mesons are heavy, the effect of binding by a 
confining potential is not negligible at presently accessible values of the photon 
virtuality Q 2 . The binding effect influences very strongly the Q 2 dependence 
of the production amplitude. 

• Even for high values of Q 2 the production of transversely polarized vector 
mesons receives important contributions from regions where the two produced 
quarks are widely separated. 

• In spite of a clear transition from the perturbative to the nonperturbative 
regime, there are nevertheless striking systematic features in the production of 
light and heavy vector mesons, whose understanding requires the incorporation 
of nonperturbative methods. 

• The gluon distribution of the proton as obtained from DIS allows only to cal- 
culate production amplitudes with zero momentum transfer. This requires 



quite essential extrapolations in the analysis of the experimental data. Since 
the approach discussed here is based on a space-time picture it takes skewing 
effects automatically into account, and therefore the dependence of the pro- 
duction amplitude on (moderate) values of t, the momentum transfer from the 
virtual photon to the proton, can be calculated. 

• In QCD the production of vector mesons is closely related to purely hadronic 
scattering processes without requiring the use of vector dominance models. 
In order to obtain a unified picture, a handling of nonperturbative effects is 
clearly necessary. 

• Closely related to the previous item is the relation between QCD and Regge 
theory. In order to investigate possible bridges between the two approaches 
again nonperturbative methods must be used. 



The calculation of photo and electroproduction presented in this paper uses a general 
method for high energy scattering based on the functional integral approach to QCD 
and on the WKB method, which is capable of incorporating both the perturbative 



and nonperturbative aspects 







121 ] . The nonperturbative input is given by a special 
model of nonperturbative QCD, called stochastic vacuum model jlH ll^ . that has 
been successfully applied in many fields, from hadron spectroscopy to high energy 
scattering. The very satisfactory results of the model in purely nonperturbative 
regions gives a strong weight to our determination of nonperturbative correction 
terms near the perturbative regime. The en ergy dependence is based on the two- 
pomeron model of Donnachie and Landshoff |l5| . 

In previous papers we have investigated photoproduction and electroproduction 
|l7| of vector mesons. We have also used the same framework to investigate 
some general features of photo and electroproduction of all S-wave vector mesons 
[rsj l which arise from the structure of the overlaps of photon and vector meson wave 
functions. 

In the recent years many more data have been obtained in HERA experiments, with 
higher accuracy and statistics, and their comparison with theoretical calculations 
provide opportunities to understand and describe important general features of the 
dynamics governing these processes. In particular, we may learn in what extent 
the experimentally observed features are contained in the global nonperturbative 



aspects of the systems in the final and initial states, such as their wave functions 
and the long range correlation properties of the intervening forces. In this paper 
we describe very successfully most of these recent data, using the same framework 
that has been tested in several other cases, without the introduction of any free 
parameter. 

Our paper is organized as follows. 

In Sec. 2 we describe the methods used in the theoretical calculations of photo and 
electroproduction of vector mesons. Since this has already been done on several oc- 
casions, we only give a short and schematic description, and then focus the attention 
on the comparison with the usual perturbative treatment. 

In Sec. 3 we present the results of our calculations of photo and electroproduction of 
all lS-wave vector mesons and compare them with the experimental data. We show 
that the predictions for all observables, the absolute values as well as the dependence 
on the photon virtuality Q 2 , the momentum transfer t and the energy W are very 
satisfactory. 

In the summary we comment on our results and present them in context. 

The appendix gives formulas and details of the theoretical calculations and presents 
some of their general properties. 

sectionGeneral formula? 

and the strictly perturbative limit 

We start from a general approach to scattering based on functional integrals and the 
WKB approximation ^} , which has been adapted to hadron hadron scatte ring 
H Q and to photo and electroproduction of vector mesons hadrons 
The basic features can be seen in Fig. Q that represents the loop-loop scattering 
amplitude in real space time. The space time trajectory of the photon is represented 
by a quark-antiquark loop, that of the proton by a quark-diquark loop. The tran- 
sition to observable electroproduction amplitudes of hadrons is achieved through 
a superposition of the loop-loop amplitudes with the light cone wave functions of 
the hadrons and the photon used as weights. This leads to an electroproduction 
amplitude of the form 
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FIG. 1: Loop-loop scattering. 
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Here S(b, W, Zi, R±, 1/2, R 2 ) is the scattering amplitude of two dipoles with separa- 
tion vectors Ri, R 2 , colliding with impact parameter vector b; q is the momentum 
transfer vector of the reaction, with 



t = -g 2 - mW + My)/W 4 + 0(W 



-q 2 



In these expressions Q 2 = —p 2 * is the photon virtuality, W is the center of mass 
energy of the 7*-proton system, and t = (py — p 7 *) 2 is the invariant momentum 
transfer from the virtual photon to the produced vector meson; z\ is the longitudinal 
momentum fraction of the quark in the virtual photon and in the vector meson, and 
we call z\ = 1 — Z\\ tf)(z, R) represent the light cone wave functions. 

The differential cross section is given by 

1 



da 

d\t\ ~ 16vrH/ 4 ' 
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For the special case of forward production (q = 0) this approach reduces to the 



dipole model for electroproduction 



221. 



Details of the evaluation of the loop-loop amplitude in the stochastic vacuum model 
can be found in previous publications |2J| . 



A. The perturbative limits 
and nonperturbative corrections 



For the photon wave functions ipj*\(zi, R±, Q 2 ) we use the well known lowest order 
expressions, keeping the same notation used before [3, [3 • 

If we concentrate on high values of virtuality Q 2 the two valence quarks stay close 
together and we may use the dipole cross section of a small object of size R\, that 
is 

J(q = 0,W,z 1 ,R 1 ) = C(W,Q 2 )Rl . (1.6) 
If we furthermore ignore the transverse extension of the vector mesons, replacing 

ip V \(zi, Ri) -> tpvxizu 0), (1.7) 
we can perform the R\ integration in Eq. ()1.2|) explicitly. 

Later we introduce specific ansatze for the vector-meson wave functions ipv\{ z i, Ri) 
and use the results of a specific model, the stochastic vacuum model, for the evalu- 
ation of the loop-loop scattering amplitude J(q, W, Z \, R\), but for the moment we 
stick to the general formulae and first study the corrections to the strictly pertur- 
bative limit. 

We first consider the simpler case of longitudinally polarized photons, where we 
obtain with the replacement (jl.7j) (for later discussion we keep the dependence on 
the quark mass m/ in our expressions): 



\/Vv f 1 8-7T 7 2 7 2 

T^ Vr , x = =-(2iW*)l 6 C(W,Q*)e f ^Q I { —^^ ^ u 0) . 

(1.8) 

The wave function at the origin ipy (zi, 0) is related to the coupling of the vector 
meson to the electromagnetic current fy by 

fv = e v V3—^= f dz x 162i%0yo(2i,O) . (1.9) 

V47T J 

Here e/ is the quark charge, and ey is the effective charge in the meson, that is 
ey = l\/2 for the p and ey = l/3\/2 for the cu meson , while for the 0, ip and T 
mesons we have ey = if. 

We introduce 

Vl(Q 2 ) = (QV4 + m)f f Q d Zl (^t+^2) 2 V [ dz x z&^zu 0) (1.10) 
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and then write the perturbative expression as 

T^t Vp>x =o = (-4^) f v ^Vl(Q 2 ) {q2/a Q + m})2 C(W, Q 2 ) . (1.11) 

Noting that the longitudinal wave function is suppressed at the end points z\ = 0, 
z\ = 1, we see that in the limit Q 2 — > oo the mass rrif can be neglected against 
ziZiQ 2 and the expression r]L becomes independent of Q 2 

«i«iVvo(«i,0) 

This is the correction factor due to the longitudinal extension of the meson as 
obta,ned in Q (up to a facto, 2 doe to different definitions). However, fo r practical 
purpose we keep using the expression (jl.K)j) with finite quark masses. If there were 
no binding effects, Z\ would be 1/2 and therefore r] L would approach 1 for Q 2 going 
to infinity. 

Comparing Eq. (jl.llj) with the perturbative expression for the longitudinal electro- 
production of vector mesons (^J we obtain 

C(W, Q 2 ) = ^-xG(x, Q 2 )a s {Q 2 ) , (1.13) 

where x is the Bjorken variable x = (Q 2 + M v )/W 2 . This yields the well known re- 
lation between the dipole cross section and the gluon density j^. After introducing 
our specific model for the meson wave functions in subsection II CI we will present nu- 
merical values for ^(Q 2 ), which represents the longitudinal momentum correction 
to the pure perturbative calculation of the amplitude for longitudinal photons. 

Nonperturbative effects also lead to a finite extension of the vector meson and 
to a deviation of the simple quadratic behaviour of the "dipole cross section" 
J(q = 0, W,Zi,Ri) in Eq. (jl.fi)) . We present numerical values for the rather large 
corrections due to these effects in subsection II CI 

As it is well known, the treatment of the transverse polarisation is more delicate 
for at least two reasons. 

• There is no strong suppression of the photon wave function at the end points 
Z\ — 0, Z\ — 1 and therefore the effective scale, namely Z\Z\ Q 2 might be 
quite low even for highly virtual photons. This makes among other things 
a systematic 1/Q 2 analysis impossible since the factor corresponding to Eq. 
(jl,12j) diverges due to the singularities at the end points Z\ — 0, 1. 



• The meson wave function is supposed to be more complicated and relativistic 
contributions are likely to be important even for rather heavy mesons. 

For simplicity we assume that the wave function of the vector meson has the same 
tensor structure as the transverse photon. Then the overlap function (jl.2j) brings 
with the replacement (jl.7|) the form 

/fa f 1 , 4mJVVi(*i,0) + 16 to 2 (z 2 + z?)V$(*i,0) 



- 2iW 2 C(W,Q 2 )e f ^P f d Zl 

Z7T JO 



2tt Jo "~ 1 {z^Q 2 + m}) 2 

where 

The relation to the decay constant fy is here given by 

, . V6 1 f\, m 2 4 !vi (z u 0) + 2uj 2 (z 2 + z 2 )^l(z u 0) 

f y = ev T7 — 7j= / dzi = • ( L16 

My a/47t Jo Z\Z\ 

This leads to a relativistic correction factor t]t which is of order 0(u 2 /m 2 ), given 
by 

(n^ (n*iA^ 2,2 f 1 , 4m| ^(^,0) + 16 u 2 (z 2 + g)j>gj(^0) 

i 4m^i(^ 1 ,0) + 8^ 2 ( 2; 2 + ^)^ 1 J( 2;i ,0) 

d^i — * — = . (1.17) 

o 4^i2i 

The final result for the for transverse polarisation is written analogous to Eq. , 
that is 

T^ty^ = (-4.Z W 2 ) f v ^VTm {Qy4 Q +mj)2 C(W, Q 2 ) . (1.18) 

Only in the absence of binding effects it would be u = 0, Z\ = ~, and we would 
then have ^r(Q 2 ) — ► 1 for Q 2 — > oo, as it is in the case of longitudinal polarization. 
Numerical values for the correction coefficients tjt, Vl evaluated in a specific model 
are given in subsection II CI Strong deviations of t]l and 7]? from unity indicate 
important deviations of the amplitudes from the pure perturbative calculation, due 
to longitudinal momentum in the photon-meson overlap and to relativistic effects. 

Important additional corrections arise from the finite extension of the vector mesons, 
and they are also discussed in subsection II CI 
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B. The specific non-perturbative model 



In order to obtain numerical information on the corrections discussed in the two 
previous subsections, we need to use models both for the wave functions and for the 
interaction of gluons and hadrons. 

The photon enters through its usual perturbative wave function. In model calcula- 



tions 



241 ] it has been shown that the perturbative wave function can be used also for 



small values of Q 2 + m 2 if an appropriate constituent mass is introduced. For light 
mesons, the masses have been determined by comparison with the phenomenolog- 
ical two point function for the vector current. The approach has been successfully 
applied in the theoretical calculation of structure functions, Compton amplitudes 
and photon-photon scattering. 

uding rela- 
2J. As 



21 



The wave function has been adapted from photons to vector mesons inc 
tivistic corrections motivated by the structure of the vector current Il7 . 

nn 

in previous papers jlfiL I17I| . we use two forms of meson wave functions: the Bauer- 
Stech-Wirbel (BSW) [25! ] . and the Brodsky and Lepage (BL) j3] forms, which 
are detailed in the Appendix. We use this type of wave function for all vector 
mesons, with quark masses determined from a best fit to the vector current (they 
are m u = = 0.2, m s = 0.3, m c = 1.25, m;, = 4.2 GeV) . 

Our nonperturbative treatment of the high energy process is based on functional 



integration 
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been described in several occasions 



13 



14|. The method has 



and on the stochastic vacuum model 

27 . 2s| . and we onh 



characteristic features. The model is based on the assumption that nonperturbative 
QCD can be approximated by a Gaussian process in the coulour field strengths; the 
gluon field correlator is therefore the quantity determining the full dynamics. Its 
parameters are taken from lattice calculations (2^. The model yields confinement 
in non-Abeliangauge theories and leads to realistic quark- ant iquark potentials for 
heavy quarks [2a] ■ It can be used to determine the loop-loop scattering amplitudes 
mentioned in the introduction and visualized in Fig. ^ 

For the energy dependence we have introduced in the model [3, 0, 0] the two- 
pomeron scheme of Donnachie and Landshoff [3] . Small dipoles couple to the hard 
and large dipoles to the soft pomeron. The transition radius was determined through 
the investigation of the proton structure function. Again we refer to the literature 
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for more information and collect some details and the relevant parameters in the 
Appendix. 

C. Numerical results 

for the nonperturbative corrections 

In order to exhibit the importance of the nonperturbative contributions we dis- 
play their effects explicitly in this subsection. We hasten to add that for the final 
calculations as given in Sec. |H] we do not split our results into perturbative and 
nonperturbative parts, but give directly the full theoretical results. 

The correction factors t]l (11.10)1 and tjt (11.17)1 are displayed in Fig. [2] as functions 
of Q 2 + 4m 2 for several vector mesons. The factor r]L reflects the effect of the 
distribution of the longitudinal momentum; we notice that it remains of order 1, and 
its influence is rather weak. In contrast, the correction in the transverse amplitude 
due to the factor tjt, reflecting mainly relativistic corrections to the transverse wave 
function, is very important, especially for high values of Q 2 + m 2 . Its large values at 
high Q 2 for the production of p-mesons indicate that the measurement of the ratio 
of the longitudinal to the transverse production cross sections tests mainly the wave 
function. According to our calculations, only the quantity tjt for light mesons is 
very sensitive to the specific choice of the wave function. 

There are large effects due to the finite extensions of the mesons. We denote by E the 
ratio of the full amplitudes (|1.2)1 to the purely perturbative ones (ll.ll)) . ()1.18)) .The 
ratio E is displayed in Fig. El left, the results for both polarisations are very similar. 

The effects are by no means negligible, even for the heavy mesons. The suppression 
due to the finite extension at low values of Q 2 leads, in the full range of presently 
available data, to a much weaker decrease in Q 2 than inferred from purely pertur- 
bative expressions. 

For large values of the quark-antiquark separation Ri, the dipole cross section differs 
from the pure i^-behaviour as given by (fb6|l . (fTT8)l . We denote by D the ratio of 
production amplitudes of the more realistic dipole cross section obtained with the 
stochastic vacuum model divided by the result obtained with the purely quadratic 
expression (jl.6j) . This factor D is also displayed in Fig. El right, we see that it is 
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FIG. 2: Correction factors r/^ Ql.lUJl and t]t, (j!.17() due to the longitudinal momentum 
distribution in the longitudinal and transverse wave function of the mesons. The factors 
are similar for the BL and BSW wave functions, except for t}t for light mesons. The 
differences in these cases are illustrated by the comparison of the long-dashed and short- 
dashed curves corresponding to the BL (long-dashed) and the BSW (short-dashed) wave 
functions for the p-meson. 
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FIG. 3: Correction factors in the evaluation of amplitudes for electroproduction. The 
factor E in the left is due to the finite extension of the meson wave function and the 
correction factor D in the right is due to nonperturbative corrections to the simple i?f 
dependence of the dipole cross section. The factors are very similar for longitudinal and 
transverse polarisations. 

important only for the light mesons at low values of Q 2 . Also D is similar for the 
longitudinal and transverse polarisation. 

The importance of nonperturbative contributions to vector meson photoproduction 



has been stressed before 
tions. 



32|], using a different approach to nonperturbative correc- 
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II. THEORETICAL RESULTS 

AND COMPARISON WITH EXPERIMENT 



A detailed description of the framework of our calculations can be found in our 
previous papers llll 13 I2I I . The basis for the determination of the loop-loop 
scattering amplitude represented in Fig. [H is the model of the stochastic vacuum. 
The calculations are determined by two parameters, the correlation length of the 
nonperturbative gluon fluctuations and the strength of the gluon condensate, which 
can be extracted from lattice calculations 

All parameters have been determined from other sources, and they are collected 
in Appendix A. Our calculations contain therefore no adjustable parameter and all 
theoretical results are true predictions. 

In the following we present a comparison of our theoretical results with the available 
data of elastic photo and electroproduction of the S-wave vector mesons, namely 
J/i/;,r,p,u,(j). 

A. Photo and electroproduction of the J/ip meson 

In Fig. we show the data for the integrated elastic cross section a of J/ip elec- 
troproduction at the fixed ene rgy W=90 GeV as a function of the photon virtuality 
Q 2 . The data are from Zeus Q| and HI j^l collaborations. The theoretical 
calculations (solid line) are made using the BL (Brodsky-Lepage) form of the vec- 
tor meson wave function; the BSW-wave function yields very similar results. The 
agreement between data and the theoretical calculations is remarkable. The dotted 
line in the figure, which nearly coincides with our theoretical calculation is a fit to 
the data proposed in the experimental paper, with the usual form 

.4 



a 



;i + Q 2 /M^) n 



(2.1) 



Fig. El shows the ratio R of the longitudinal to the transverse cross sections 



R = ^f, (2.2) 
a 1 



again for J/ ip elastic electroproduction at W=90 GeV. As the ratio cancels influences 
of the specific dynamical model, this quantity tests directly details of the overlaps 

12 



10 



10 



— 

c 



10 



y* p -» J/\)/ p 




■ HI 2005 
• Zeus 2004 
o Zeus 2002 

solid: theoretical calculation 



dotted: fitting by experimentalists 
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FIG. 4: Integrated elastic cross section of J/ip photo and electroproduction at the en 



34( and HI 



ergy W=90 GeV as a function of Q 2 . The data are from the Zeus 
collaborations. The solid line represents our theoretical calculation using the BL wave 
function. The dotted curve is a pure fit to the experimental points of the form given by 
Eq.flHH). 



of wave functions, helping t 
The data are from Zeus 134 



ongitudinal and transverse structures 
]| collaborations at HERA. As seen 



le study of their 
"1 and HI Q 

in the figure, the presently available data are not very accurate. The theoretical 
calculations with the two kinds of wave function - BL and BSW - give nearly the 
same results for a = a L + a T and for R = a L /a T . 

Fig. El shows the comparison of our calculations for the energy dependence of cross 



sections with the recent Zeus and HI data |3 
data from fixed target experiments 38, 



34 |35| and older photoproduction 



at lower energies. For illustration of our 
results for electroproduction we show data and theoretical results in the Q 2 range 
6.8 - 7 GeV 2 for which there are both Zeus and HI data. 

Often the W dependence of experimental cross sections is fitted through single 
powers in the form 

a = Const, x W 5{q2) , (2.3) 

which may be useful in limited W ranges. Values of the parameter 5 obtained 
at several values of Q 2 in J/ip electroproduction are compared to our theoretical 
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FIG. 5: Ratio R of longitudinal and transverse cross sections for J/ 



34 



electroproduction 



3d and HI 



35 



as function of Q 2 for fixed energy W=90 GeV. Data from Zeus 
collaborations at HERA. The solid line and dashed lines show our theoretical calculations 
respectively with the BL (solid) and the BSW (dashed) wave functions. A good numerical 
representation for the BL result is R{Q 2 ) « 0.76 (Q 2 /M 2 )/(l + Q 2 /M 2 ) om . 

predictions in Fig. [?l based on the two-pomeron scheme. 

In contrast to purely perturbative approaches, our theoretical treatment allows to 
calculate the dependence of the cross section on the momentum transfer t. The 



35| at 



data on the t-distribution in J/ ip photo |33|, |35[ and electroproduction 
W=90 GeV are shown in Fig. |H] together with our theoretical results (solid lines). 
Also shown are electroproduction data and calculations at Q 2 = 6.8 — 7 GeV 2 where 
measurements from both Zeus and HI are available. 

Our calculation predicts a curvature in the log plot for the t distribution, which we 
may describe with the form 



da 
d\t\ 



da 



d\t\ 



xF(\t\) 



t=o 



da 
d\t\ 



-b\t\ 



x 



t=o 



[i + a\t\y 



(2.4) 



For 7*p — > pip at W=90 GeV our calculations at Q 2 = give a = 4.06 GeV -2 
and b = 1.75 GeV -2 . Present data give no clear cut evidence for this curvature. 
The distribution becomes flatter as Q 2 increases, but in the investigated range does 
practically not change with the energy (no shrinking). 

In order to determine the integrated cross section the experimental data have to be 
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FIG. 6: Energy dependence of the integrated elastic cross section for J ftp production 



from the HERA collaborations 



33 



34 



35j |. compared with our theoretical calculations. 



38, 



39]. 



The fixed target photoproduction data are from the E-401 and E-516 experiments 
Electroproduction is represented by the Q range 6.8-7.0 GeV 2 , for which there are both 
Zeus and HI data points. 



extrapolated to the point of minimum momentum transfer. This brings a theoretical 
bias to the experimental points. It is therefore meaningful to compare our theoretical 
results with observed data at small but finite momentum transfer. This is done in 
Fig. El the agreement is excellent. 



B. Photo and electroproduction of T meson 



For the T there are only two data points for photoproduction, at < W >=120 



GeV (3] and < W >=143 GeV 41 1, shown together with our theoretical energy 
dependence in Fig. I n the second plot we show our calculation of the Q 2 

dependence at the fixed energy W=130 GeV, together with the data for Q 2 = at 
the energies 120 and 143 GeV. Both plots show a fairly good agreement with the 
data. It should be noted that, in spite of the high mass scale, the finite extension 
of the meson has a considerable effect, as can be seen in Fig. Elleft. 
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FIG. 8: The t dependence of the differential cross sections of J/i/j elastic photo 



and electroproduction 
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35| |. The Zeus(2002) photoproduction results have separate 
information from events with ip decaying into (full circles) and e + e~ (empty circles). 

The solid lines represent our theoretical results. 
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FIG. 9: Q 2 and W dependences of the (nearly) forward differential cross sections of 

\ fin 

35l | and electroproduction [341 . 1351 ] . The solid lines represent our 



J/ip elastic photo 
theoretical results. 
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FIG. 10: Integrated elastic cross section of T photoproduction as function of the en 
ergy and electroproduction at the energy W=130 GeV as a function of Q 2 . The data 
at the energies < W > = 120 and 143 GeV are respectively from Zeus 40J and HI 
collaborations. The solid lines represent our theoretical calculations using the BL wave 
function. 
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C. Photo and electroproduction of p meson 



The classical fixed target experiments 42, |43j, |44j, |45|, |46j provide important reference 
for the magnitudes of cross sections in p photo and electroproduction at center of 
mass energies about 20 GeV. The data are shown in Fig. ^2 together with the 
results of our calculation. The agreement is satisfactory given the quite important 
discrepancies within the data in the interval Q 2 = 4 — 7 GeV 2 . 



Y* p -> p p 

W=20 GeV 
fixed target experiments 
solid line: our calculation 
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FIG. 11: Integrated elastic cross section for p electroproduction at the energy W=20 GeV 
as a function of the photon virtuality Q 2 . The data are from fixed target experiments 
42, 



43 



44 



45 



461 ] . spanning almost two decades. The solid line represents the results of 



our calculations. 



Fig. H2l shows the data for the integrated elastic cross section a of p electroproduction 
at the fixed energy W=90 GeV as a function of the photon virtuality Q 2 . The data 



are from Zeus 



361 WR and HI 



48. 



collaborations. The solid line represents 
the theoretical calculations using the BL form of the vector meson wave function. 
Our theoretical calculations give an good overall description of the data but do not 
reproduce the details of the Q 2 dependence of p electroproduction as well as they 
did in J/ip production. The present data indicate that in the p meson case a fitting 
of a simple form like Eq. (j2.1j) is not satisfactory for the whole Q 2 range. Two 
of such forms may be applied, with a transition in the values of both parameters 
(normalization and power) occurring somewhere in the Q 2 range from 5 to 7. More 
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experimental measurements are needed to clarify the structure of the data in this 
region, which may contain important information about the dynamics of the process. 
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FIG. 12: Integrated elastic cross section of p elec tro-p roduction at the energy W=90 GeV 
as a function of Q 2 . The data are from Zeus jsfil. L?! and HI 4^ 4jJ] collaborations. The 



solid line shows our theoretical results using the BL wave function. 

Fig. El shows the ratio R = a L ja T of cross sections for longitudinal and transverse 
polarisations, for p elastic electroproduction at W=20 and 90 GeV. The data at 



90 GeV are from NMC 



46j , Zeus 



M |47| and HI Q, Q, Q- The w=2 ° GeV 
data are from the E-665 experiment |45j . To indicate the differences, the theoretical 
results for both BL (solid line) and BS W (dashed line) wave functions are displayed. 
For large Q 2 , R becomes very sensitive to the values of the small transverse cross 
section. The measurement of this quantity R therefore provides important tests on 
the structure of the meson wave function. 

Fig. El shows the energy dependence of p electroproduction for several values of 
Q 2 together with our theoretical results. On the right-hand-side plot we show the 
effective power 5 for the energy dependence (see eq. (j2.3j) ) for the IV-region of 
approximately 30 to 130 GeV. The experimental points are from Zeus 
and HI 0, 3- O ur theoretical description, solid line, is very satisfactory within 
the experimental errors. 

There are no published measurements of the differential cross section da/d\t\ in p 
production for nonzero values of Q 2 . The Zeus photoproduction data at W=75 and 
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FIG. 13: Ratio R between longitudinal and transverse cross sections for p electroproduc 



tion as function of Q for fixed energies W=20 GeV , with data from E-665 [45r j. and W=90 



GeV with data from NMC 



461 ] . Zeus 



36, 



4J and HI 



48 



49 



45], 



501 ] . The solid and dashed 



lines show the theoretical calculations with the BL and BSW wave functions respectively. 



94 GeV 



51 



52f are shown in Figs. E3 In the low t range there are new preliminary 



data from HI |53| at several energies. The numbers for the differential cross sections 
at 70 GeV for low t extracted from their plots are included in the figure, and they 
seem to confirm the previous Zeus measurements at W=75 GeV 5Jj. Our theoretical 
calculations are shown in the figures, the description is satisfactory. 



D. Photo and electroproduction of to meson 

As particles of about the same size and mass, u and p have similar behaviour in the 
soft processes that we study here. 

The Q 2 dependence of u photo and electroproduction 54 , |2] is shown in Fig. 



The agreement between our results and the data is not perfect, but satisfactory, in 
view that there is no free parameter involved in the calculations. 



The existing data 



56 



571 ] on the t dependence of the differential cross section in 
uj production are shown in Fig. El I n the figure are put together the data points 
of the energies W=15 GeV (with Q 2 = 0) and W=80 GeV (with Q 2 =0.1 GeV 2 ), 
and the corresponding theoretical curves. The coincidence of shapes exhibits the 
universality of the form factors of t dependence in our model. Our prediction of a 
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FIG. 14: Energy dependence of p electroproduction cross section in the region from 
approximately 30 to 130 GeV. The solid lines show our theoretical calculations for several 
values of Q 2 and the data are from Zeus 



36, 



42L 



51( and HI 



48 



collaborations. The 



plot in the right hand side shows data and theoretical values for the parameter 5 of the 
energy dependence W s . 



curvature in the plot of da/d\t\ seems to be confirmed by the data. 



;he integrated elastic cross section in 6 
S0fi, » with our theoretical 



E. Photo and electroproduction of meson 

The <j) meson has a strategic place between the J/ ip and the p meson and may help 
to understand the differences in behaviour of heavy and light vector mesons and 
also to clarify the interplay of perturbative and nonperturbative aspects of QCD. 

In Fig. El we show the Q 2 dependence of 
photo and electroproduction at W=75 GeV 
calculations. As in the case of p electroproduction, there is indication that the data 
cannot be well represented by a single expression of the form of Eq. (|2.1j) . The data 
are poorer here than in the p case, and it is important to investigate the possibility 
of a transition region in an intermediate Q 2 range below 10 GeV 2 in which the 
normalization and the power change values rather rapidly. 

The data for the ratio R = a L /a T for W=75 GeV as function of Q 2 j^, 0] are 
shown in Fig. EI together with our results with BL and BSW wave functions. The 
calculations exhibit, as in the p meson case, the sensitivity of the ratio R to details 
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FIG. 15: Experimental data and our theoretical prediction (solid line) for the t dependence 
of p photoproduction cross sections. The published data are from the Zeus collaboration, 
at the energies 75 GeV (low t) 



51| and 94 GeV (large t) 



£2j. The 94 GeV data are 



0.16 



rescaled with a factor (90/94) 
preliminary information from HI 



0.993 in the figure. We also include in the low-t figure 



53j, extracted from their plots, for W=70 GeV. 
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FIG. 16: Q 2 dependence of the cross section of oj elastic production 
line represents our calculations. 



0,Q 



. The solid 



of the wave functions. Offering a reference for the two kinds of calculation, the plot 
also shows (dotted line) a fit of the data, made by experimentalists. 

The effective power 5(Q 2 ) describing the energy dependence as W s in <p electro- 
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FIG. 17: t dependence of a; photoproduction cross sections |56l. 1571] . The plot puts together 
data points at the energies W=15 GeV (with Q 2 = 0) and W=80 GeV (with Q 2 =0.1 
GeV 2 ), and the corresponding theoretical curves. The curvature and the similarities of 
shapes of the t distributions for two different energies and Q 2 values are characteristic 
features of our framework. 

production j^] is shown in Fig. EH an d compared with our results using the two- 
pomeron model. The theoretical Q 2 dependence is similar to that of J/ip and p 
production, whereas the experimental points, with large errors, indicate a flatter 
behaviour. Obviously more data are necessary. 

Our calculations give very good descriptions for the Q 2 dependence of the differential 
elastic cross section in forward directions in <j) electroproduction, as it does in the 
J/ip case. Fig. |2~T1 shows the comparison with the experimental data [sol ] . 

There are no published measurements of da/d\t\ in <ft production for nonzero Q 2 . 



The Zeus photoproduction data at 94 GeV |K 



are shown in Fig. 1221 Our 



theoretical calculations give a very satisfactory description. 



F. Results concerning several vector mesons 

The quantitative predictions made in a unique way for different kinds of vector 
mesons cover several scales of magnitudes in cross sections. This global coverage is 
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FIG. 18: Q 2 dependence of (j) production cross sections at W=75 GeV and our theoretical 
description. The data are from Zeus 
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60j and HI 



59]. 
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FIG. 19: Q 2 dependence of the ratio R = a L /a T in (f> electroproduction at fixed en- 

ergy W=75 GeV. The solid and dashed lines represent our calculations with BL (solid) and 
BSW (dashed) wave functions; the dotted line is a fit of the form R = 0.51 (Q 2 /M 2 )°- m . 



exhibited in Fig. [221 The same global description is given for the forward differential 
cross section shown in Fig. EU In this figure the charge factors squared ey, see Eq. 
(II. 9j) . for each kind of meson are extracted, making the quantities almost universal in 
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FIG. 20: The effective power 5 parameter governing the energy dependence of <f> electro- 
production. The data are from Zeus |60j, the solid line is our theoretical result. 
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FIG. 21: Q 2 dependence of the differential cross sections of <fi elastic in the forward 
direction. The solid lines represent our calculation. The second plot contains the same 
information, showing the typical behaviour of a straight line in the variable Q 2 + My. 



a Q 2 + My plot. However, the universality is only approximate, and our calculation 
predicts correctly the observed displacements. 

The data on the energy dependence of the integrated cross sections for p, <fi and i/j 
mesons have been presented and compared to the theoretical predictions for each 
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FIG. 22: t dependence of </> photoproduction cross sections and its theoretical description 
in our nonperturbative calculation with the stochastic vacuum model. The published data 



are from the Zeus collaboration, at the energy 94 GeV, first at low t 



Q 



58| and then at larger 



case. The parameter S(Q 2 ) of the suggested simple energy dependence 

a(Q 2 ) = Const, x W sm 

has also been given in each case. This parametrization is an approximation valid 
in a limited energy range, since the true energy dependence in our calculation is 
determined by the two-pomeron scheme, but it is considered useful in practice. 

We then evaluate 5 using the energy range W=20 - 100 GeV, for all values of Q 2 . The 
results are put together in Fig. We note that all curves start at the minimum 
value 4 x 0.08 at the same unphysical point Q 2 + M v = and all are asymptotic to 
5 = 4 x 0.42 as Q 2 increases. We then have the form of parametrization 

(1 + Q 2 /M 2 ) n 



0.32 + 1.36 



A + (1 + Qi/Ml 



(2.5) 



and the values for A and n are given in table □ 



Our model has definite predictions for the t dependence in differential cross sections. 
The shape can be conveniently represented by the form given in Eq. (|2.4jl . Fig. l2*6l 
shows the form factor .F(|£|) for all vector mesons at fixed energy W=90 GeV. 
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FIG. 23: Integrated elastic cross sections for all vector mesons at W=90 GeV, as functions 
of Q 2 + My. The lines represent our theoretical calculations. 
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FIG. 24: Differential cross sections in a forward direction for all vector mesons at W=90 
GeV, as functions of Q 2 + My , with extraction of charge factors . The lines represent the 
theoretical calculations with the stochastic vacuum model as in Fig. I2.il 



The curvature in the log graph is a prediction of our framework. The values of 
the parameters as functions of Q 2 are shown in one of the plots. In Table |H] some 
numerical values of a, b are given. 
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FIG. 25: Parameter 5(Q 2 ) of the energy dependence of the cross sections. 

As the virtuality Q 2 grows, the ranges of the overlap functions decrease, and the 
electroproduction cross sections of all mesons become flatter, all tending together to 
the shape characteristic of the T meson, with same limiting values a = 4.02 GeV 2 
and b = 1.60 GeV 2 for the parameters. The limiting shape of the distribution for 
very large Q 2 is illustrated in Fig. |2H where we see all vector mesons superposed. 
In the figure we draw the bit of straigth line representing the slope considered as 
the average for the interval from \t\ = to \t\ = 0.2 GeV 2 . As indicated inside the 
plots, the calculations of form factors presented in the figures are made for W = 90 
GeV. In the second plot presented in Fig. EH we show the (absence of) dependence 
of the form factor on the energy W. Thus we predict that there is no shrinking of 
the forward peak in da/d\t\ as the energy increases. The experimental data are not 
yet sufficient to test all these predictions. 
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FIG. 26: Form factor of the \t\ distribution in the elastic differential cross sections 

photoproduction of vector mesons, with the characteristic curvatures in the log plot pre- 
dicted in our calculations. The curves follow the shapes given by Ea. (|2.4|) . with a(Q 2 ) 
and b(Q 2 ) represented in the second figure. Specific numerical values of the parameters 
for Q 2 = , 10 and 20 GeV 2 are given in Table ITT1 The values in the limit of very large 
Q 2 are a = 4.02 GeV -2 and b = 1.60 GeV -2 . 



TABLE II: Some values of the parameters a and b of Jfe 
The full Q 2 dependence is shown in Fig. 1261 
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p(770) 
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4.349 
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4.145 


1.846 


cj(782) 


7.20 


3.08 


4.323 


1.990 


4.141 


1.835 


</>(1020) 


5.40 


2.77 


4.211 


1.934 


4.091 


1.807 


J/^O-S) 


4.06 


1.75 


4.026 


1.696 


4.022 


1.670 


T(1S) 


4.03 


1.61 


4.027 


1.606 


4.025 
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III. SUMMARY AND DISCUSSION 



We have shown the predictions for elastic electroproduction processes using two ba- 
sic ingredients: 

1) the overlaps of photon and meson wave functions, written as packets of quark- 
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FIG. 27: The plot in the left hand side shows the shape of the \t\ distribution in the 
differential cross section for electroproduction common to all vector mesons for very large 
Q 2 . The straight line passes through the points \t\ equal to zero and 0.2 GeV 2 , indicating 
the average slope that would be measured in this limit. The plot in the right hand side 
shows that i*X|i|) does not depend significantly on the energy, and thus measurement of 
slopes at fixed Q 2 would give a constant value for each vector meson. 

antiquark dipoles, with protons described also as packets of dipoles (in a convenient 
diquark model for the nucleon), 

2) the interaction of two dipoles described in terms of geometric variables in an 
impact parameter representation of the amplitudes based on nonperturbative prop- 
erties of the QCD gluon field. 

These quantities put together and integrated over the distribution of dipoles in 
initial and final states lead to a correct description of the data concerning all vector 
mesons. 

In all cases the variations with energy are very well described by the Regge picture, 
with soft and hard pomerons coupled to large and small dipoles, respectively. We 
recall that in approaches mainly based on perturbation theory the energy (or x) 
dependence is introduced through the gluon distribution in the proton. 

Each of the different mesons enter in the calculation characterized only by the masses 
and charges of its quark contents, and with their normalized wave function individ- 
ualized only by the corresponding electromagnetic decay rate (related to the value 
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of the wave function at the origin). 

The specific nonperturbative input is the stochastic vacuum model, which has been 
successfully applied in many fields, from hadron spectroscopy to high energy scat- 
tering. The basic interaction of two dipoles depends only on universal features of 
the QCD field, which are the numerical values of the gluon condensate and of the 
correlation length of the finite range correlations. These two quantities have been 
determined by lattice investigations and tested independently in several instances 
of phenomenological use of the dipole-dipole interaction. 

As has already been pointed out the main features of the Q 2 dependence of 
electroproduction of vector mesons are contained in the overlap integral. For values 
of Q 2 attainable at present this overlap integral is determined by perturbative QCD, 
through the photon wave function, and by nonperturbative QCD, through the me- 
son wave function. The production of transversely polarized mesons is determined 
by the meson wave function even for very high values of Q 2 . The importance of 
nonperturbative effects in vector meson production has also been stressed in j^, 
where the perturbative two gluon exchange has been supplemented by an exchange 
of nonperturbative gluons. In our approach, which incorporates both perturbative 
and nonperturbative effects, we are able to give a fair overall description of all ob- 
servables of vector meson production: the energy dependence, the Q 2 dependence, 
the ratio of longitudinal to transverse mesons and the angular distribution. 

In this paper our calculations are compared to the large amount of recent HERA 
data for p, u, 0, J/0 and T mesons. 



APPENDIX: QCD AND WAVE FUNCTION PARAMETERS 



Here we recall some expressions and the numerical values of quantities used in the 
nonperturbative contributions, which have been unchanged for many applications. 

The loop-loop scattering amplitude S(b, W, z\, R±, z 2 = 1/2, R 2 ) determining the 
essential quantity J(q, W, z\, R±) in Eq. (J1.3J) . can be calculated using the stochas- 
tic vacuum model. The essential input parameters of this model are the correla- 
tion length a of the gauge invariant two gluon correlator and the gluon condensate 
(g 2 FF). These quantities have been determined in lattice calculations (3). The 



31 



numerical values used in this and previous papers are 

a = 0.346 fm {g 2 s FF) a 4 = 23.5 . (A.l) 

The constant k = 0.74 appearing in the correlation functions preserves its value 
determined by lattice calculations. 

The energy dependence is based on a two-pomeron model, small dipoles with size 
R < R c couple to the hard pomeron with an intercept aph(0) = 1.42, whereas large 
dipoles with R > R c couple to the soft pomeron with an intercept ap s (0) = 1.08. 
The transition radius R c has been determined from the ^-dependence of the proton 
structure function to R c = 0.22 fm. 

For the proton wave function occurring in Eq. (jl.3|) we use a diquark-quark Gaussian 
wave function with the transverse radius Rp = 0.75 fm. 

The light quark masses which can simulate confinement effects in the otherwise 



perturbative photon wave function 



2J| are 



m u = md = 0.2 GeV m s = 0.3 GeV . (A.2) 



For the heavy quarks we take the renormalized masses in the MS scheme 

m c = 1.25 GeV m b = 4.2 GeV (A.3) 

The full form of photon and vector meson wave functions used in our work, including 
the helicity dependences, have been presented before jl^ . Q] , with two forms for the 



vector mesons: the Bauer-Stech-Wirbel (BSW) 25j, and the Brodsky and Lepage 
(BL) 2(| wave functions, where the separate r, z dependences are respectively 



V 



>BSW{z,r) = —j== Vzz exp 

V47T 



Mi, 1. 



2u 2K 2 

(here My represents the vector meson mass) and 



exp[-^ 2 r 2 ] , (A.4) 



>BL{z,r) = —==exp 

V47T 



m)(z-\) 2 



exp [-2 zzu 2 r 2 ] , (A.5) 



2zzu 2 

with rrif representing the quark mass. 

The parameters N (normalization) and u> (that fixes the extension) are determined 
using the electromagnetic decay rates of the vector mesons. Their values are col- 
lected UM in Table EH 
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TABLE III: Parameters of the vector meson wave functions 









BSW 






BL 








transverse 


longitudinal 


transverse 


longitudinal 


Meson 


JV (<^eV) 
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0.2618 


2.0469 


0.3088 


1.8605 


0(1020) 


0.07592 


0.2568 


4.6315 


0.3549 


4.6153 


0.3113 


1.9189 


0.3642 


1.9201 


J/^(1S) 


0.26714 


0.5770 


3.1574 


0.6759 


5.1395 


0.6299 


1.4599 


0.6980 


2.3002 


T(15) 


0.23607 


1.2850 


2.4821 


1.3582 


5.9416 


1.3250 


1.1781 


1.3742 


2.7779 



After summation over helicity indices, the overlaps of the photon and vector meson 
wave functions 

PrvAQ 2 ; zi, R-0 = Vva^i, Ri)*^ r a (<5 2 ; «i, Ri) (A. 6) 

that appear in Eq. ([1.2j) are given by 



Prv,±i;BSw{Q 2 ;z,r) =e v — — (e f u 2 r[z 2 + z 2 \K x {e f r) +m 2 f K (e f r)) <f) B sw{z,r) 

Z7T v ' 

(A.7) 

and 



P-y*v,±i;BL{Q 2 ]z,r) =e v -^-(4ef uj 2 rzz[z 2 + z^K^ej r)+m 2 f K (e f r)) (f) BL (z,r) 

(A.8) 

for the transverse case, BSW and BL wave functions respectively. For the longitu- 
dinal case we can write jointly 

\/3a 

Prv,o-,x(Q 2 ;z,r) = -16e v — — u z 2 z 2 Q K (e f r)<f> x {z,r) , (A. 9) 
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with X standing for BSW or BL. The effective quark charges ey are l/y/2 for p, 
1/3 \/2 for lu, —1/3 for and T and 2/3 for ip. The quantity 



e f = yjz(l -z)Q 2 + m 2 (A.10) 
and the modified Bessel functions are introduced by the photon wave functions. 
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